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ABSTRACT 

In this work, we identify carbodiimide (HNCNH), which is an isomer of the well- 
known interstellar species cyanamide (NH2CN), in weak maser emission, using data 
from the GBT PRIMOS survey toward Sgr B2(N). All spectral lines observed are in 
emission and have energy levels in excess of 170 K, indicating that the molecule likely 
resides in relatively hot gas that characterizes the denser regions of this star forming 
region. The anticipated abundance of this molecule from ice mantle experiments is 
~10% of the abundance of NH2CN, which in Sgr B2(N) corresponds to ~ 2 x 10 13 
cm -2 . Such an abundance results in transition intensities well below the detection limit 
of any current astronomical facility and, as such, HNCNH could only be detected by 
those transitions which are amplified by masing. 
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1. Introduction 

Historically, searches for new astronomical molecules resulted in the detection of favorable, 
high line strength transitions based on a thermal approximation to the excitation of these species 
in interstellar environments. At the temperatures of hot molecular cores inside molecular clouds, 
these transitions often reside at (sub)millimeter wavelengths. However, line confusion near the 
Boltzmann peak can lead to ambiguous identifications, and the peak intensities of low abundance, 
large organic molecules may never rise above the noise floor or the confusion limit. The success of 
molecule searches at centimeter wavelengths has shown that unique excitation conditions can lead 
to the unambiguous identification of very low abundance species with high accuracy. 

Carbodiimide (HNCNH) is the seco nd most stable isomer of cyanamide (NH2CN ), with cyanamide 



being more stable by ~4.0 kcal mol 1 ( Duvernav et al. 20041 ; iTordini et al. 20031 ). Since the de 



tection of NH2CN by Turner et al. (1975) toward the high-mass star-forming region Sgr B2(N), 
HNCNH has been recognized as a candidate interstellar molecule. HNCNH is also a tautomer of 
NH2CN, and can be formed via an isomerization reaction whereby a H atom can migrate along the 
molecular backbone from the amine group. At room temperature, HNCNH is in thermally-induced 
equilibrium with NH2CN at a fraction of Given typical abundances of NH2CN and tem- 

peratures in many astrophysical environments of << 300 K, gas-phase tautomerization of NH2CN 
is unlikely to produce HNCNH in appreciable abundance. Yet, HNCNH may exist in detectable 
abundance out of thermal equilibrium with NH2CN as in the case of interstellar hydrogen cyanide 
(HCN) and interstellar hydrogen isocyanide (HNC) - in which HNC is in greater abundance in 
some clouds (see Allen, Goddard, & Schaefer 1980 and references therein). 

Although tautomerization is likely inefficient in the gas phase under interstellar conditions, ex- 
perimental studies have shown NH2CN — > HNCNH conversion in water ices and matrices is far more 
efficient (Duvernay et al. 2005 and references therein). In the solid phase, the association of up to 
five water molecules with NH2CN on an ice surface has been shown to sign ificantly lower the act i- 



vation barrier to tautomerization and promote the production of HNCNH ([Duvernav et al. 20051 ). 
In fact, HNCNH formation on a water-ice surface is shown to occur at temperatures as low as 70 
K. The relative abundance of HNCNH formed from thi s process has been me asured to range from 



4% of NH 2 CN at 80 K to as much as 13% at 140 K (jDuvernay et al. 20051 ). In the laboratory 



H NCNH sublimation occ urs between 80 - 170 K, with no additional desorption observed above 170 



K (jDuvernay et al. 20041 ). In the ISM, non-thermal desorption via shocks is also a likely liberation 
mechanism. Assuming that tautomerization of NH2CN on dust grain ice mantles is the dominant 
formation pathway for HNCNH, and that subsequent desorption occurs for both species, the gas 
phase abundance of HNCNH may be capped at ~10% of NH2CN. 



NH2CN has an estimated column density of ~2x 10 14 cm -2 towards Sgr B2(N) ( Nummelin et al. 2000i ): 
we therefore anticipate an HNCNH column density on the order of (~10 13 cm -2 ). If thermal emis- 
sion at hot core temperatures dominates the spectrum of HNCNH, the most favorable high line 
strength transitions are at millimeter and submillimeter wavelengths. Nevertheless, we have pre- 
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viously used the Robert C. Byr d Green Bank Te lescope (GBT) to detect n ew molecular species 
including trans- methyl formate (jNeill et al. 20121 ) . and cyanoformaldehyde ([Remijan et al. 2008l ) 
which only had measurable astronomical line intensities detected at centimeter wavelengths. En- 
couraged by these results, we searched for centimeter wave transitions of HNCNH. All data were 
taken as part of the PRebiotic Interstellar MOlecular Survey (PRIMOS), and targeted rotation- 
torsion transitions of HNCNH selected from the Cologne Database for Molecular Spectroscopj0EI 
(jMiiller et al. 20051 ). Our results are presented in §2. In §3, we evaluate the probability of mis- 
assigning these features to HNCNH and discuss evidence that observed transitions are masing. 
Further, we explore the usefulness of our technique as a method for new molecule detection. 



2. Observations and Results 

All data used for this project were taken as part of the National Radio Astronomy Observatory's 
(NRAO) 100-m Robert C. Byrd Green Bank Telescope (GBT) PRebiotic Interstellar MOlecule 
Survey (PRIMOS) Legacy Projectlf] This NRAO key project started in Jan 2008 and concluded in 
July 2011. The PRIMOS project recorded a nearly frequency continuous astronomical spectrum 
from 1 to 50 GHz towards the Sgr B2(N) molecular cloud, with the pointing position centered on 
the Large Molecule Heimat (LMH) at (J2000) a = 17h47ml9.8s, 5 = -28°22' 17". An LSR sourc e 



velocity of +64 km s" 1 was assumed. Intensities are presented on the T^ scale ( Ulich &: Haas 19761 ). 
See Neill et al. (2012) for further details of the observations, full data reduction procedure, and 
analysis. 

Targeted transitions of HNCNH are shown in Table Q] which lists each transition's rotation- 
torsion doublet transition quantum numbers, calculated rest frequency (MHz), logarithm of the 
Einstein A coefficient, lower state energy (cm -1 ), and Gaussian fitted astronomically measured 
line intensity (mK) and FWHM line widths (km s" 1 ) in the observed spectra. The table provides 
Gaussian fit parameters for peak intensity and line width for components at vlsr — +64 km s _1 and 
+82 km s . Molecular material at these two velocities are characteristic of clouds that lie within 
the GBT beam along the same line of sight (see e.g. Remijan et al. 2008b and references therein). 
Figure 2 shows the spectra for the observed passbands. The red and blue vertical lines indicate the 
location of the emission feature at LSR velocities of +64 km s _1 and +82 km s _1 , respectively. We 
find four unblended lines of HNCNH in this region (panels a) - c) and g)). One frequency range was 
unobservable by the GBT (16 GHz) while two spectral feature were contaminated by transitions 
of CH3OH (panels d) & g)). Finally, one transition was not observed due to lack of frequency 
coverage (46.2 GHz). Of the observable transitions, we show clear detections of both the +64 km 



1 Original laboratory data cataloged from Birk et al. (1980), Wagener et al. (1995), & Jabs et al. (1997). 
2 Available at www.splatalogue.net ( Remiian et al. 2008h . 



3 Access to the entire PRIMOS dataset and specifics on the observing strategy including the overall frequency 
coverage, is available at |http://www.cv.nrao. edu/~aremijan/PRIMOS71 
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s and +82 km s components at 4.3, 4.8 and 25.8 GHz, and the +64 km s component at 45.8 
GHz (Figure [I]) . The four line detections and two non-detections are consistent with only masing 
lines being detectable, as discussed in §3. 
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Fig. 1. — Carbodiimide (HNCNH) spectral passbands toward Sgr B2(N) recorded from the GBT 
PRIMOS Survey. Rotation-torsion doublet transition quantum numbers are shown in each panel. 
The passband width displayed is 500 km s~ x in each case. The spectra are plotted as a function 
of frequency (MHz), corrected for a LSR source velocity of +64 km s . The blue and red vertical 
lines indicate the location of the transition rest frequency (see Table 1) at an assumed LSR source 
velocity of +64 km s~ 1 and +82 km s -1 , respectively. Data in all panels were Hanning smoothed 
for display purposes. 
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Table 1. HNCNH Rotation- Torsion Doublets, Spectroscopic, and Observed Astronomical 

Parameters 
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a All lines except those at 4 GHz have been experimentally measured with (Obs - Calc) frequency uncer- 
tainties < 30 kHz. U ncertainties on the 4 GHz lines are calculated to be < 30 kHz, type A, k = 2 (2a) 
dTavlor fc Kuvatt 1994h - 



b The uncertainties for the intensities and linewidths are type B, k = 1 (lc) ( Taylor Sz Kuvatt 1994 )- 
c No data 

d Transition contaminated by CH3OH emission. 

f Estimated based on a gaussian fit to the low-frequency shoulder of the CH3OH emission at +82 km s _1 
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Fig. 2. — Energy level structure of the relevant transitions for the 4.8 GHz maser (top) and 36.6 
GHz line (bottom). Energy levels are ordered by increasing energy, but are not drawn to scale. 
Allowed transitions are indicated by arrows. 
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Table 2. Spectroscopic Line Parameters for the Transitions Labeled in Figure [2] 
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3. Discussion 

The initial identification of HNCNH was based on only the two lines at 4.3 and 4.8 GHz. In 
order to quantify the probability of coincidental overlap of features at these frequencies which could 
lead to a possible mis-assignment, we used the method outlined in Neill et al. (2012). We find 
37 observed transitions within a representative window of 200 MHz of PRIMOS at C-band (4-6 
GHz), 14 of which are within ±50% of the intensity of the detected features. This line density is 
typical of the PRIMOS survey in this frequency range. Of these 14 features, 4 were in emission. 
If we then assume a conservative FWHM line width of 25 km s^ 1 (almost twice our measured 
FWHM) we calculate the probability of a single line falling coincidentally within one FWHM of 
our line to be 0.75%. For two detected lines, that probability drops to 0.002%. We find this is both 
compelling evidence for a clear detection and an excellent example of the power of the GBT in new 
molecule detections at centimeter wavelengths, where line confusion is drastically reduced relative 
to the millimeter and sub-millimeter regimes. 

Initially, the large beamwidths and high-energy nature of the observed lines pointed to emission 
arising from the hot, extended source that surrounds the cold LMH and couples well to the beam at 
these frequencies. An LTE analysis, however, showed excessively high column densities for thermal 
emission arising from these lines. In order to further verify the identification, and to understand the 
unusually high intensities of the lines, we examined higher frequency transitions. Unfortunately, 
the GBT has no receiver that operates between 15-16 GHz; we detected an emission feature at 
25.0 GHz, but its doublet transition at 25.5 GHz is blended with the 9(2,7)-9(l,8) transition of 
CH3OH at 25.541 GHz; we detected no emission features at 36.6 and 37.1 GHz; and a feature was 
identified at 45.8 GHz, the highest frequency line currently covered in the PRIMOS survey. The 
lack of observed transitions at 36.6 and 37.1 GHz was initially troubling as there has been a firmly 
established series of requirements for an interstellar detection, including that no transition which 
should be of an observable intensity be "missing" (Snyder et al. 2005). However, an examination 
of the rotational energy level structure of HNCNH offers an explanation. 

Figure [2] shows energy level diagrams containing transitions that are relevant to the observed 
masing lines at 4.8 GHz (labeled "M", Figure 2(a)) and 36.7 GHz (labeled "X", Figure 2(b)). The 
quantum mechanical parameters associated with transitions in Figure 2 are given in Table El In 
the case of the 4.8 GHz line "M", the upper energy level (A2 S 18(0,18)) is coupled to higher energy 
levels via fast millimeter and sub-millimeter transitions; the only downward transition out of this 
level is very slow, with an Einstein A coefficient of ~ 10~ 9 s" 1 . Downward and upward transitions 
out of the lower level (Aid 17(1,17)), labeled "A - D", however, occur rapidly with Einstein A 
values between 0.2 - 2 x 10~ 2 s" 1 . This creates a population inversion between the A2 S 18(0,18) 
and Aid 17(1,17) levels, leading to the observed masing at 4.8 GHz. An analogous process occurs 
for the torsion doublet, causing masing at 4.3 GHz. A similar energy structure exists for the 25.5 
GHz and 46.7 GHz transitions and their torsion doublets, leading to masing in these transitions as 
well. 
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For the 36.7 GHz line (labelled "X", Figure 2(b)), the upper state (A 2d 15(1,15)) is again 
coupled to higher energy levels via fast transitions, with the 36.7 GHz transition being the slow 
emission path to A 2s 16(0,16). In this case, however, a fast downward transition is allowed to the 
A 2s 14(0,14) state (transition labelled "N"). Population inversion is not achieved because of the 
presence of this drain, and the 36.7 GHz transition does not mase. A similar scenario holds for the 
15.9 GHz transitions, as well as their torsion doublets. In summary, an analysis of the energy levels 
of HNCNH predicts a population inversion resulting in masing for the torsion doublets at 4, 25, 
and 46 GHz, with no masing by the doublets at 15 and 36 GHz. Furthermore, the 4 GHz doublets 
should show a greater degree of inversion due to the longer lifetime of the transitions (Table [T]) 
than the 25 or 46 GHz transitions, resulting in more intense lines, which is in excellent agreement 
with our observations. 

The observed maser emission is relatively weak and displays no line narrowing, suggesting 
that these are unsaturated maser lines. As collisional and radiative pumping rates for HNCNH are 
unknown, a full pumping model was not possible and is beyond the scope of this work. A zeroth- 
order calculation of the critical density for the fast transitions connecting the maser transitions 
was, however, performed. Using collisional coefficients for similar transitions in HNCO, we find 
critical densities required to thermalize these transitions on the order of 10 8 - 10 10 cm~ 3 , far higher 
than those typically associated with the region. Given this, and that the Sgr B2(N) region is 
known to contain a broad range of favorable radiative excitation conditions, a radiative excitation 
mechanism for these transitions is possible. The current observations, however, do not provide 
a definitive answer to the excitation question, and further investigation is warranted. A mapping 
study of the 4 GHz lines in Sgr B2(N) would be extremely beneficial to elucidating the environment 
these signals are arising from, including their spatial correlation with other known maser such as 
CH3OH, and therefore the most likely excitation mechanisms. 

We have also carried out an analysis of the expected intensity of HNCNH lines in this frequency 
region based on a thermal population of HNCNH. Assuming local thermodynamic equilibrium 
(LTE) conditions, we calculate a column density of HNCNH giving rise to the 4 GHz maser lines 
of ~2 x 10 16 cm" 2 at 150 K, although this value is largely invariant over temperatures from 80 
- 500 K. Given this column density, the expected intensity of the 36.7 GHz transition would be 
740 mK. Further, at a column density of ~2 x 10 13 cm" 2 (~10% of NH2CN in this source 



(jNummelin et al. 2000 )). we would expect the most intense lines in this region to be less than 1 
mK for rotational temperatures above 80 K. Given this, and the lack of any detection at 36.7 GHz 
with an RMS noise level of ~10 mK, we conclude these transitions are not arising from a thermal 
population. 

This is consistent with our non-detection of transitions which are not masing. In gas above 
80 K, the most intense lines of HNCNH at LTE would fall at millimeter and sub-millimeter wave- 
lengths. Numerous high-resolution, high-sensitivity spectral line surveys of this region have been 
carried out (see Wang et al. 2011, Tercero et al. 2010 & Nummelin et al. 2000 and refs. therein) 
covering the these transitions up to ~1.8 THz. With an abundance of 10% NH2CN, the strongest 
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transitions under LTE emission at 80 K would be on the order of the RMS noise level of the survey 
observations. With increasing temperature, the expected line strengths decrease rapidly. Thus, at 
these column densities, a thermal population of HNCNH would likely not be detectable in even the 
most sensitive line surveys to date. 

This is in itself remarkable. A heroic effort has been performed in recent years on compiling 
molecular inventories and facilitating new molecule detections using state-of-the-art millimeter and 
sub-millimeter observatories and molecular line surveys. These surveys have provided invaluable 
information on the physical and chemical conditions present in the ISM and their high scientific 
value is unquestionable. They suffer, however, from extremely high line densities and long integra- 
tion times, relative to lower-frequency telescopes. Because of this, identifications of new molecular 
species, especially those in low abundance, can be difficult due to the line confusion at the noise 
floor and the degree of coincidental overlap of target lines with other molecular transitions. 

In this work, we have shown that a molecule which would be undetectable via transitions arising 
from a thermal population can be identified via maser transitions at centimeter wavelengths. This 
detection was possible only through low frequency surveys of a chemically rich region by the GBT 
through a dedicated project (PRIMOS). The number of transitions which display significant maser 
activity in any given molecule is likely to be s mall, but weak masing behavior has been observed in a 



number of molecules such as H2CO and NH3 (jForster et al. 1980l : lGaume et al. 199ll ). In addition 



at centimeter wavelengths the lack of line confusion makes definitive identification of a species 
possible with only a small number of observed transitions. This may represent a new strategy for 
searches for key molecules of interest to the astro chemistry and astrobiology communities which 
have not yet been detected due to their low abundances. 

In summary, we have detected four transitions of carbodiimide (HNCNH) towards Sgr B2(N) 
with very high confidence. All four signals have been found to be the result of maser activity. We 
also report two transitions that were not detected, consistent with HNCNH only being detectable 
towards Sgr B2(N) through maser lines. This detection presents a new methodology for searches 
for interstellar molecular candidates which may be too low in abundance to be detected in thermal 
emission by modern radio observatories. 
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